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Abstract 

Giant phospholipid vesicles are shown to exhibit thermally acti- 
vated transitions between a prolate and an oblate shape on a time scale 
of several seconds. From the fluctuating contour of such a vesicle we 
extract ellipticity as an effective reaction coordinate whose temporal 
probability distribution is bimodal. We then reconstruct the effective 
potential from which we derive an activation energy of the order of 
ksT in agreement with theoretical calculations. The dynamics of this 
transition is well described within a Kramers model of overdamped 
diffusion in a bistable potential. Thus, this system can serve as a 
model for macroscopic bistability. 

PACS: 05.40+j, 68.15+e, 82.70-y. 
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Macroscopic bistable systems for which thermal energy is sufficient to 
overcome the barrier in observable times are rare, quite in contrast to mi- 
croscopic bistability so ubiquitous in nature |]T|. In this paper, we consider 
bistability in the shape space of giant (lOfim) phospholipid vesicles 0. The 
shapes of these vesicles are determined by the minima of the bending en- 
ergy of the vesicle membrane subject to constraints on the total area and 
enclosed volume || U A cut of the energy landscape along a path con- 
necting nearby local minima then defines an effective double well potential 
in shape space. Since this potential is created by the vesicle due to its bend- 
ing elasticity, our system is qualitatively different from the most prominent 
example of macroscopic bistability studied so far which is a one micron silica 
bead caught in a double well potential generated by an optical trap 0. In 
the latter system, the (optical) potential is created externally in real space 
via gradients of the electromagnetic field. 

A vesicle (shape) is not a static entity but shows thermally excited fluc- 
tuations, which are clearly visible by video microscopy, due to the softness of 
the vesicle membrane. At constant external parameters such as temperature, 
these fluctuations occur around a well-defined mean shape corresponding to 
the global minimum of bending energy. As temperature changes, these mean 
shapes usually evolve smoothly except at phase boundaries denoting so-called 
discontinuous transitions along which shapes of different symmetry have the 
same energy. In this report, we discuss the observation of vesicles at the 
prolate-oblate transition |7], g]. These vesicles exhibit thermally activated 
jumps between the two states at constant external parameters. In Fig. 1, we 
show the contours of a typical sequence of vesicle shapes at equidistant times 
(6.3 s). A vesicle is seen to change its shape from a prolate with the mean 
symmetry axis within the focal plane to an oblate with the mean symmetry 
axis perpendicular to the focal plane and vice versa [[|. 

For a quantitative analysis, the contours of the vesicles in the focal plane 
are determined by image analysis [IUL ITTL |T2|. We describe the shapes using 
polar coordinates, 

r (4>) = r \l + J2 a n cos(n0) + b n sin(n0) , (1) 

\ n>2 n>2 J 

where ai = h\ = implies the definition of the center of gravity of each 
contour as the origin. Since gravity holds the vesicles on the bottom of the 
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Figure 1: Prolate-oblate fluctuations of a DMPC vesicle. The snapshots 
are taken at equidistant time intervals of At = 6.3 s. The vesicle is seen 
to change its shape from a prolate with the mean symmetry axis within 
the focal plane to an oblate with the mean symmetry axis perpendicular to 
the focal plane and vice versa. Vesicles were prepared from stearoyl-oleoyl- 
phosphatidylcholine (SOPC) or dimyristoyl-phosphatidylcholine (DMPC) in 
sucrose solution (50 mOsm) and incubated in excess glucose solution (48 
mOsm). They were observed by video phase-contrast microscopy in a mi- 
crochamber at constant temperature. Due to the slightly higher density of 
sucrose, the vesicle sink to the bottom of the observation chamber and rest 
gently against the glass wall. Thus, the free rotational diffusion of vesicle 
shapes is restricted to the 2 dimensions of the focal plane of the microscope. 
This stabilisation procedure is crucial for recording time series of vesicle 
shapes. Suitable vesicles were selected and observed near the boundary be- 
tween the prolate and the oblate phase. 
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Figure 2: Time series of vesicle ellipticity A = ya^ + bl near the prolate 
oblate transition. The ellipticity exhibits fluctuations around two distinct 
mean values which correspond to the prolate and oblate shape respectively. 
The shapes are quantified by extracting the vesicle contour via image anal- 
ysis. Briefly, a contour point is defined by the gradient of the phase con- 
trast profile across the membrane with respect to the local background level. 
Thereby, one finds a closed loop made up of points with sub-pixel resolution. 
The amplitudes a2 and 62 are then obtained from this contour via Eq. ([!]). 



chamber, oblate shapes have a circular contour while prolate shapes (with 
their major axes parallel to the bottom) look ellipsoidal |^3j. Thus, it is 



sufficient to monitor the amplitude A = J (a^ + > 0> which effectively 
measures the ellipticity of the vesicle contours. 

In Fig. a time series of A for a typical vesicle is shown. The apparent 
jumps in A between a small numerical value near zero and a larger positive 
number reflect the bistable nature of the system, already evident in Fig. |l[ 
The larger values correspond to the fluctuating prolate shape, whereas the 
small values parallel the oblate vesicle with its fluctuating equatorial circular 
contour. 

The diffusion in shape space is evident in Fig. |3|, where we show the 
distribution of (02,62)- The points cluster at the center (the oblate phase) 
and along a circle (the prolates) because of the rotational symmetry of the 
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Figure 3: Distribution of instantaneous contour parameters (02, 62) i n a 
polar diagram. The dark center corresponds to the oblate shape, the darker 
rim to the prolate phase. Since the activation energy is of the order of k B T, 
there is significant population also along the saddle points. 
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Figure 4: Histogram of events A of the ellipticity obtained with a bin size 
0.045. The dotted lines are guides to the eye. The double well potential 
V(A) is then reconstructed by using Eq. fl2|). The saddle point location at 
A s = 0.21 and its curvature wl = 29 fcgT are obtained via a fit (solid line) 
of the saddle to a parabola. 

problem. 

We now construct an effective potential for the variable A from the mea- 
sured probability distribution P(A), represented by a histogram of A as 
shown in Fig. [|. The probability to find the shape in the interval [A, A + dA] 
(and any orientation a) is proportional to dA Ae~ v( - A ^ kBT ', where V(A) is an 
effective potential. The additional factor A arises from changing from carte- 
sian coordinates a 2 , 62 to polar coordinates A and a. From the measured 
P{A), we thus get 

V(A) /k B T = In A - In P{A) + const (2) 
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as shown in Fig. The activation energies Q b,pro from the oblate and the 
prolate phase, respectively, thus obtained, are Q a b — 1.5 ±0.1 ksT and 
Q pro ~ 0.2 ±0.1 k B T. 

Although the prolate shape is metastable with a small barrier height 
Qpro, h is nevertheless populated with a higher probability than the oblate 
minimum (compare Fig. 0). This is due to the measure factor Ad A for 
the probability distribution P(A), which reflects the larger density of states 
(with energy V(A)) of the prolate compared to the oblate shapes. All prolate 
vesicles with a fixed ellipticity A but varying orientation a have the same 
energy V(A), while the oblates lack this orientational degree of freedom. 

This observation of a vesicle as fluctuating in a bistable potential |14 



is qualitatively different from the three other fluctuation effects observed in 
vesicles so far. (i) Fluctuations of quasispherical vesicles have long been used 
to extract the bending rigidity of bilayer membranes |TD|, |TT|, [T5|, [TJJ. Our 
experiment shows that this technique must fail if the available excess area 
of the vesicles is so large that the discontinuous character of the prolate- 
oblate transition becomes relevant, (ii) The budding transition |TJ|, [T7|, 
where a small satellite is expelled from a parent vesicle, is a strong first-order 
transition for which the activation energy cannot be surmounted thermally. 
One then observes pronounced fluctuations as the spinodal is approached 
[ O] . (iii) Conformal diffusion of vesicles of higher topology involves diffusion 
along a one parameter continuous curve in a "flat" potential |19|, ^0 . 



A quantitative comparison of our experimental results with theoretical 
first-principle calculations is delicate. First, the presence of gravity makes 
even a calculation of the shape of lowest energy non-trivial [fL3|l . However, 
since we expect gravity to be a minor correction of the energy, one can, in 
a first approximation, compare with the energy of free vesicles. Ignoring all 
shape fluctuations, this energy has been calculated for the two local minima 
(prolate and oblate) as well as for the intermediate stationary saddle point 
which is a non-axisymmetric ellipsoid [p], Pf. One finds a weak first-order 
transition with a maximal energy barrier between prolate and oblate for rea- 
sonable material parameters of the order Q° ob pro — 1 ksT. In principle, this 
energy barrier acquires an additional entropic contribution from the avail- 
able phase space for fluctuations orthogonal to the reaction co-ordinate A. 
While for axisymmetric shapes, these entropic corrections can be calculated, 
a calculation for the non-axisymmetric saddle shape is far beyond presently 
available techniques. Since the measured activation energies are in reasonable 
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agreement with the theoretical value, our results indicate that these entropic 
corrections do not differ strongly for the three stationary shapes involved. 

The dynamical aspects of this transition can be modeled as an over- 
damped Kramers escape problem in a potential V(a 2 ,b2) = V(A) J21], ^2 



The diffusion in the angle co-ordinate can be integrated out, which leads to 
an one-dimensional diffusion process for the variable A in the effective po- 
tential Vrfj(A) = V(A) — k B TlnA. The corresponding Langevin equation 
reads jd t A = — dV e g/dA + (, where ( is a stochastic force with the usual 
correlations (C(^)C(O) = ^k B T^5{t — t'). The effective friction coefficient 
7 is assumed to be independent of the shape and given by 7 = const wR 3 , 
where n is the viscosity of water and R the apparent radius of the oblate. 
The Kramers mean escape time from the oblate minimum becomes [52] 



lV27rk B T 

Tab = o T~ eX P Qob T, (3) 

where the frequencies 0J o b,s are defined by a quadratic approximation to the 

2 

energy in the oblate (V(A) « V(0) + ^fA 2 ) and at the saddle at A s (V(A) » 

V(A S )-4(A-A s f). 

The pre-exponential factor in (|3|) becomes more transparent if it is ex- 
pressed as t ob y / 2irk B T/ui s A s in terms of measurable quantities. Here, t a b = 
7/^06 is the relaxation time for a 2 or 6 2 fluctuations in the oblate phase. 
The measured dynamical correlation functions (a 2 (£)a2(0)) and (& 2 (t)fo 2 (0)) 
can both be fitted against simple relaxational behavior with a decay time 
t b — 5 ± 2 s. Using the measured values Qob/ksT ~ 1.5 , A s ~ 0.21 
and oj 2 s ~ 29 k B T as extracted from the potential (see Fig. f|), we obtain 
r b ~ 50 ± 25 s. This value is in reasonable agreement with the mean resi- 
dence time f ob ~ 76 ± 8 s in the oblate phase as extracted directly from the 
time series A(t) as shown in Fig. ^|. 

We now turn to a discussion of the prolate phase. For the fluctuations 
of the ellipticity A around its mean value (A), we find a relaxation time 
ta — 8±ls. Fluctuations in a amount to rotational diffusion. For the latter, 
we find a rotational diffusion coefficient D = (5.7±0.2) x 10~ 3 rad 2 / s. It is this 
free diffusion which is responsible for the relative stability of the metastable 
prolate state. The mean escape time from the prolate phase should not be 
written in a Kramers form since the implicit assumption of Kramers theory 
that the activation energy is larger than k B T is violated on the prolate side. 
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To summarize, we have reported the observation of a macroscopic bistable 
system as evidenced by large shape fluctuations of giant vesicles between a 
prolate and oblate state. This result confirms the theoretical prediction of 
a weak discontinuous transition between prolate and oblate shapes. The 
mean escape time out of the metastable oblate state was found to be in 
reasonable agreement with the value obtained within a Kramers approach for 
the effective radial potential. Thus, fluctuating vesicle shapes may serve as a 
model system for further studies of macroscopically slow diffusion processes 
in a bistable potential. 

We are grateful for discussions with E. Evans, M. Jaric, M. Kraus, R. 
Lipowsky, and M. Wortis. 
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